Temperature dependence of magnetism near defects in SrB 6 . 
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The T-dependence of magnetic moments in SrB6 is stud- 
ied through spin-polarized band calculations for a supercell 
of Sr27Bi56 containing a B6 vacancy. The magnetic moment 
decays rather quickly with T despite the fact that only elec- 
tronic Fermi-Dirac effects are included. This result and the 
T-dependence of moments near a La impurity can hardly ex- 
plain the reports of a very high Curie temperature in hexa- 
borides. but suggest that the magnetism is caused by some 
other type of impurity. 

Early reports of weak ferromagnetism (FM) in lightly 
doped CaBg or SrBg, associated with an unusually high 
Curie temperature (Tc) 1 , have been followed by several 
theoretical and experimental studies. For instance, it 
has been speculated that the FM could be of excitonic 
origin 2 , be a spontaneous polarization of the dilute elec- 
tron gas 3 , or be caused of special conditions around impu- 
rities and thereby lead to band magnetism 4 ' 5 . The band 
structure of the pure hexaboride 6 , which is metallic when 
calculated within the density functional (DF) approach, 
has been questioned, since one GW-calculation found a 
gap at the Fermi energy, Ep 7 . These results are in con- 
trast to a recent independent GW-calculation by Kino et 
al, where the ground state is metallic 8 . Also, there are 
conflicting results concerning the interpretation of the 
experimental results on the weak magnetism. Quantum 
oscillations indicate that the FM state is not intrinsic to 
the bulk, but that impurities might be important 9 . This 
conclusion was supported by Mori and Otani 10 who found 
that the FM signal could be removed by a treatment of 
the crystals, and suggested that that the FM is due to 
iron impurities. Other defects such as vacancies have 
been suggested as well 11 . Possible vacancies of whole 
Bg-clusters were anticipated in calculations by Monnier 
and Delley 12 who showed that a sizable FM moment sur- 
rounded such a vacancy in a supercell of CaBg. They 
suggest that the observation of weak ferromagnetism in 
dilute doped hexaborides is caused either by a surface or 
interface effect, since the Bg vacancy is modeling the nat- 
ural cleavage planes in polycrystals, or by some amount 
of Bg vacancies in pure crystals. 

In the present letter, by use of the linear Muffin-Tin 
orbital (LMTO) band method in the local spin-density 
approximation, we study the variation of the magnetic 
moment around a Bg vacancy in Srg as function of tem- 
perature. This is done by self-consistent spin-polarized 
calculations, in which T-dependent electronic excitations 
are modeled by the Fermi-Dirac distribution, but no ther- 
mal disorder of the lattice is taken into account. As in 



ref. 12 we consider a 3x3x3 supercell with one Bg vacancy 
in SrBg so that the total cell is S^Bisg. The missing 
Bg atoms form empty MT-spheres, and the calculations 
include additional empty spheres in the open part struc- 
ture. The cell contains 19 nonequivalent MT-spheres, 4 
Sr, 8 B and 7 empty ones. Other details of the calcula- 
tion are as in ref. 5 . The density-of-states (DOS) functions 
near Ep for the non-polarized case, are shown in Fig. 1 
for two sets of k-points (4 and 10) in the irreducible part 
of the Brillouin zone. It is seen that the essential fea- 
tures of the DOS are developed already when using 4 
k-points, with only small differences in the fine details. 
In particular, the DOS peak at Ep is similar in the two 
sets of k-points, so that the criterion for Stoner mag- 
netism should come out well already when using only 4 
k-points. The Stoner factor is calculated to be 1.05, just 
above the limit for ferromagnetism, 1.0, indicated by the 
horizontal line in fig. 1. The broken line shows the partial 
DOS from the six equivalent B-atoms closest to the Bg- 
vacancy. It is seen that roughly 1/3 of the total DOS at 
Ep comes from these sites although they represent only 
about 3 percent of the total number of atoms in the cell. 
In the spin-polarized calculations (at low T), when using 
4 k-points, one finds a magnetic moment, m —2.1 fj,B per 
cell, in agreement with ref. 12 . This moment decreases to 
about 1.8 /is when the number of k-points is increased 
to 10 for the same temperature. This type of magnetism 
is a standard Stoner magnetism 13 , where the large DOS 
on the 6 B-atoms nearest to the vacancy contribute most 
to the Stoner instability, cf. Fig 1. Similarly, about 40 
percent of the total moment in the spin-polarized calcula- 
tions comes from these six B-sites closest to the vacancy. 

Spin-polarized calculations at different temperatures 
in the Fermi-Dirac distribution, show that the moment 
goes to zero when T is increased to about 700 K. The two 
sets of k-points have been used at the lowest and highest 
T, and both sets give a vanishing moment at T=700K. 
Calculations for intermediate T, using 4 k-points, show 
first a gradual decrease of m at low T, while the curve 
drops sharply as T is approaching 700 K, as shown in Fig 
2. This T-dependence is understood from the fact that 
Ep falls on a narrow peak in the DOS, so that thermal 
smearing, due to electronic excitations, can reduce the ef- 
fective DOS and the Stoner factor. Since Tc in the doped 
hexaborides is of the same order, 700 K, it is tempting 
to assign the observed Tc to the mechanism of Stoner 
magnetism due to the Bg vacancy. However, the typi- 
cal behavior for ferromagnets is that the temperature at 
which S goes below the critical value for ferromagnetism, 
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is generally much larger than the Curie temperature (Tc) 
for real materials 14 . By taking into account additional 
smearing effects on the DOS coming from thermally in- 
duced disorder of the structure one obtains more realis- 
tic values of Tc in mean- field calculations 15 . The atomic 
displacements from thermal vibrations at 700 K can be of 
the order 0.25 a.u. (corresponding to 7-8 percent of the 
B-B distance) for a material with a Debye temperature of 
350 K, while large smearing effects on DOS peaks appear 
already at a lower degree of disorder 16 . The present cal- 
culations include only the T-dependence via the Fermi- 
Dirac distribution and the effects from thermal disorder 
are expected to reduce the moment further. Therefore, 
it is difficult to understand the high Tc in the real mate- 
rials from the T-dependence of the moments within the 
Stoner mechanism. One possibility would be that the lo- 
cal structure around a B@ vacancy somehow can escape 
effects of thermal disorder, in which case it is justified 
to consider only the effect Fermi-Dirac smearing. The 
position of Ep relative to the peak appears to be at the 
optimal place for the largest possible Stoner factor for 
this supercell, and an increased height of the peak at 
larger T would be opposite to the normal T-dependence 
of the DOS. 

Another mechanism leading to magnetism is due to a 
gain of Coulomb energy in addition to the exchange en- 
ergy, by means of charge transfers. This mechanism can 
be activated if there is a large derivative N' of the DOS 
of an impurity band, leading to very small moments near 
the impurity 4 . This difference with the normal Stoner 
mechanism, suggests that this weak magnetism could re- 
sist better to temperature and broadening of the DOS by 
disorder, since N' can be large within a wider energy in- 
terval compared to that of a peak in the DOS. However, 
the results from calculations for a La impurity in a 3x3x3 
cell give only a partial support to this expectation 5 . The 
DOS of this cell is such that N' is maximal a bit below 
Ep, and even if there is a small moment at 600 K, it tends 
to zero at low T. It is difficult to have a sizable moment 
both at low and high T. By optimizing the position of 
N' relative to Ep (and the amplitude of N) by a com- 
bination of the virtual crystal approximation and a non- 
conventional choice of linearization energies in the band 
calculation, it is possible to have a weak moment within 
some range of T. Such band results are not ab-initio, 
but they are useful to show a correlation of moment and 
charge transfers. They also show that magnetism within 
this mechanism depend on fine details of the electronic 
structure near the impurity, which in turn suggest that 
such magnetism should be a result of circumstances. Fur- 
thermore, it is expected that thermal disorder should be 
as important for this type of magnetism as it is for the 
standard Stoner magnetism. 

In conclusion, the band results confirm that weak mag- 
netism is possible, either as standard Stoner magnetism 
around a Bg vacancy or assisted by Coulomb energies 
near a La-impurity, but it is more difficult to understand 
the high Tc values from both of these mechanisms. These 



results and recent experimental reports 9,10 motivate sim- 
ilar electronic structure studies of the properties near iron 
impurities in these materials. However, the exact state of 
iron contamination, impurity sites or clustering on sur- 
faces, is yet not clear. 
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FIG. 1. Total DOS of Sr 27 Bi56 near E F , using 4 and 10 
k-points. The broken line is the partial DOS on B 6 closest 
to the vacancy. The vertical line indicate the limit for Stoner 
magnetism. 
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FIG. 2. Calculated magnetic moment (hb per cell) as func- 
tion of temperature for Sr27Bi56. 
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